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Ps) We investigate methods for modelling metabolism within populations of cells. Typi- 

cally one represents the interaction of a cloned population of cells with their environment 
as though it were one large cell. The question is as to whether any dynamics are lost 

^ by this assumption, and as to whether it might be more appropriate to instead model 

each cell individually. We show that it is sufficient to model at an intermediate level of 

\^ granularity, representing the population as two interacting lumps of tissue. 
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1 Introduction 



The emerging field of systems biology seeks to reconcile subcellular-level components (such 
J> as enzymatic reactions) with cellular- and organism- level behaviour (such as metabolism). 

Non-linear processes dominate these interactions; experience from other areas of science has 
taught us that mathematical models, continuously revised by new information, must be used 
to describe and interpret complex biological phenomena [H |2] . 

As systems biology grows, so we see a proliferation of mathematical models of cell metabolism 
fSJ and signalling - see the many examples at the model repositories BioModels.net [3j and 



CellML.org [3]. Given the inherent difficulties in in performing single cell experiments, one 
t> property held in common by many of these models is the assumption of "lumped dynamics" . 

^ To explain this term, consider a typical scenario in which a million S. cerevisiae are grown in 

$H a chemostat. Experiments are performed to measure average metabolite concentrations over 

^ the population of yeast cells. A mathematical model of metabolism is then built in which the 

cell has these average characteristics, but a volume equivalent to a million cells (see Fig. [T]). 

Given the identical metabolic characteristics of each clonal cell, it would seem natural to 
approximate the system by lumping the population as a single mass. Intuition would suggest 
that dynamics are unchanged but, as we shall see below, this linear, verbal reasoning approach 
is incorrect. However, we show that it is not necessary to consider each individual cell - which 
would lead to a million times as many ODEs - rather correct dynamics can be captured by 
considering two interacting lumps of cells. 



*This preprint first appeared on Nature Precedings on 1 December 2009 doi:10.1038/npre. 2009. 4033.1 
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(a) 





(b) 

Figure 1: Modelling at different scales. See Eqs. ([T])-([2]) for a mathematical representation, 
where y denotes extracellular and x intracellular concentrations. Typically, (a) one models 
the cell population as one bulked compartment; at the other end of the granularity scale, (b) 
one could consider each of the n cells individually, which would lead to approximately n times 
as many ODEs. 
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2 A theorem 



We frame the problem mathematically. Let Xi denote a set of metabolite concentrations 
within cell i, and y a set of external concentrations (see Fig. [T]). Assuming each cell has 
identical characteristics, we may write 



x'i = f{xi,y) 



i = 1, . . . ,n 



y' = g{y) — ^K^i^y) 



(1) 
(2) 



Here / denotes intracellular reactions, h transport into cells and g the rate of metabolite 
supply. 

Linearise about a steady-state Xi = x*, y = y* to give stability matrix 

Ar,= : ■. : (3) 



We propose that 
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where A denotes the spectrum. That is, the system bulked into two compartments has the 
same eigenvalues as the system with three compartments, but more than the system with one 
compartment. 

To show X{An) C X(An+i), let Vn = {xi, . . . ,Xn\y)' and suppose AnVn = Xvn- Taking 



^1, . . . , Xn, / 

n ^ 



Vn^\ = \ X 

we find = Af„+i. 

Now suppose tin+i = (xi, . . . ^Xn+\\y)' and suppose = \un+\- Taking 



(5) 



nxi + Xn+1 nXn + Xn+l 



n + 1 



n + 1 



y 



(6) 



we find AnUn = An„,. 



Finally, we must consider the possibility that Un = 0, i.e. Xi = —Xn+i/n Vi. If n > 2, this may 
be overcome by first creating a new eigenvector . i = (x„+i,X2, . . . ,Xn-,xi\y) by swapping 



3 



two elements, then constructing Un as above. Thus we may conclude A(A„) ^ A(y4„+i) for 
n > 2 as required. 

The practical implication of the above theorem is that, the dynamic behaviour (or at least 
the linear dynamic behaviour) of a full system of cells may be captured by bulking the cells 
into two compartments. If cells are instead bulked as one, some behaviour will be lost. 

Moving back to specifics, we may construct the two sets of eigenvectors associated with the 
system. If ui = {x\yy is an eigenvector of Ai, then m„ = (x, . . . ,x\yy is the corresponding 
eigenvector of An- If f = x is an eigenvector of fx, then = (x, 0, . . . , 0, — x, 0, . . . , 0|0) are 
the corresponding eigenvectors of An. 



3 An example 

Prom a stability perspective, the system 

x=fix,y*) (7) 

may be naturally unstable at x = x*, but this instability may be masked in the model through 
tight control in y - leading to the eigenvalues of Ai all having negative real part. However, 
if the cells are not bulked as one, but rather as two (or more) compartments, the feedback 
exposes the realities of the system as An now inherits positive real part eigenvalues from fx- 

For example, the Brusselator is a model proposed in 1968 for an autocatalytic, oscillating 
chemical reaction [5]. In dimensionless form, dynamics may be written as 



u' = 1 - {b + l)u + au^v (8) 
v' = bu — au^v (9) 

Its steady-state is given by (n, f) = (l,&/a) and if 6 > a + 1 there exists a globally-stable 
limit-cycle (see Fig. [2]). 

This model may be transformed by setting x = (n, v) and letting y = b now be a variable 
representing the externally-supplied nutrient (similar results may be obtained by setting y = 
a). 



u'i = 1 - (6 + l)ui + au]vi (10) 

v'l = bui — au^Vi (11) 

b' = g--y\ {hiui + h2Vi + hsb) (12) 
n ^—^ 

For certain parameter values, control on b will seem to stabilise the system (n=l). (see 
Fig. [3] (a)). However, when the bulked cells are split, the underlying oscillations return (b). 
Similar dynamics are observed when comparing n = 2 and n = 3 (c). 
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Figure 2: (From Eqs. ([8])-(|9])). Stable limit cycle of the Brusselator. Parameter values used 
are a = 1, 6 = 3, u(0) = 1.01 and v{0) = b/a. 

4 Discussion 

Returning to Fig. [T| we see the two scales of granularity typically used in metabolic modelling. 
Typically one represents a population of cells as a single compartment, rather than considering 
the dynamics of n individual cells. The reasons for this are not clear. It may be that it is 
assumed that a population of clonal cells would behave in the same way as this. Alternatively, 
it may be assumed that in order to capture the interactive dynamics, around n times as many 
differential equations would be required. 

As we have shown, both mathematically and via the example of the Brusselator, neither of 
these assumptions are true. Rather, to answer the titular question, two lumps are required. 
It is hoped that by using this methodology as standard, new dynamics may be exposed that 
were previously hidden by the standard assumptions. 

Acknowledgements I acknowledge the support of the BBSRC/EPSRC Grant BB/ C008219/1 
"The Manchester Centre for Integrative Systems Biology (MCISB)". Thanks to Dave Broom- 
head for fruitful discussions. 
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(b) 





(c) 



10 20 30 40 50 60 



Figure 3: (From Eqs. ([T0])-([T2])). 
used are as in Fig. [2| with g = 2, 



(a) n = 1: steady state stabilisation. Parameter values 
hi = -4, /i2 = 0, /i3 = 2 and 5(0) = 3. (b) n = 2: stable 
limit cycle obtained by dividing populations. Parameter values used are as before, with initial 
conditions ni(0) = 1.01, ^2(0) = 0.99. (c) n = 3: initial conditions tii(O) = 1.01, n2(0) = 1 
and usiO) = 0.99. 
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